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We study the behavior of systems in which the interaction contains a long-range component that 
does not dominate the critical behavior. Such a component is exemplified by the van der Waals 
force between molecules in a simple liquid- vapor system. In the context of the mean spherical model 
with periodic boundary conditions we are able to identify, for temperatures close above T c , finite- 
size contributions due to the subleading term in the interaction that are dominant in this region 
decaying algebraically as a function of L. This mechanism goes beyond the standard formulation of 
the finite-size scaling but is to be expected in real physical systems. We also discuss other ways in 
which critical point behavior is modified that are of relevance for analysis of Monte Carlo simulations 
of such systems. 



I. INTRODUCTION 

An item of conventional wisdom in the study of critical phenomena is that the critical point behavior, including 
finite size scaling, is controlled by a relatively small number of features of a system, among which are the structure 
of the order parameter, the nature of boundary conditions, and the general properties of the interaction coupling 
the order parameter at different locations. In particular, it is believed that short-range interactions lead to universal 
critical phenomena in the case of a given system. For a non-critical 0(n) system with periodic boundary conditions, 
finite size corrections are expected to be exponentially small in the ratio L/£, where L is the smallest of the system's 
linear dimensions, and £ is the correlation length. This expectation holds everywhere on the phase diagram, with 
the possible exception of the coexistence curve, where, in certain systems, gapless spin wave excitations give rise to 
long-range correlations. 

When interactions are long ranged, the above expectation is subject to revision. The hallmark of a long-range 
interaction in the context of critical behavior is a diverging n th moment. That is, if V(r) is long-range, then, for some 
sufficiently high n, the integral 

V n = { V{r)r n d d r (1.1) 



diverges. This diverging moment appears in the Fourier transform of the interaction, v(q) through an anomaly in its 
expansion as a power series in q. In the case of a very short range interaction, the power series expansion is entirely 
in integer powers of q 2 . Any deviation from such an expansion represents an anomaly. 

That long-range interactions can alter the scaling behavior of a critical system has been known for some time . 
The first investigation of this phenomenon in the context of the renormalization group jjj established that when 
v(q) cx q a for sufficiently small q, with a < 2, then the critical point behavior of an interacting spin system differs 
fundamentally from that of a system in which the interaction is short-ranged. The upper critical dimension d u for 
any such 0(n) system turns out to be d u = 2a. This has been established by renormalization group arguments in 
and rigorously proven in ||. On the other hand, the lower critical dimension is di = a |^,||. 

In the context of critical phenomena, the criterion for short range interactions, with respect to the leading critical 
behavior, is a finite second moment of V(r). In terms of the power series expansion of v(q), this means that whatever 
anomaly exists, does not interfere with, or dominate at small q, the first two terms in the expansion in powers of q 2 . 
That is, one can write for small q 

v(q) = v + v 2 q 2 + R(q) (1.2) 



where R(q) is asymptotically smaller than the first two terms on the right hand side of ( |1.2| ) for small q. When the 
interaction between the order parameter at different points in the system has a Fourier transform that looks like 



the right hand side of (1.2), one expects that the thermodynamic critical behavior will be as predicted by standard 
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approaches for a system whose interaction is entirely short-rang ed [§. For finite -size 0(n) systems with 2 < d < 4 
this implies |^-[ic|l, among the other things, that the only relevant variable on which the properties of the finite 
system depend in the neighborhood of its bulk critical temperature T c is L/£. For T > T c . finite-size corrections for 
systems having periodic boundary conditions arc expected to be exponentially small in terms of L/£. One assumes 
that for periodic boundary conditions all reference lengths, aside from the bulk correlation length £, will lead only to 
corrections to the above picture. 

The above expectation has been challenged in recent papers by Chen and Dohm ]Tl] |l5|] . In these papers, it is pointed 
out that a model having both short-range interactions and a sharp, wavelength-dependent cutoff of fluctuations A, 
will exhibit finite size corrections to the infinite system limit that swamp those traditionally associated with finite 
size scaling. They demonstrate that for such a particular choice of the fluctuation cutoff function, one will observe 
corrections to the infinite system thermodynamic behavior going as an inverse power law in L that do depend also on 
LA and not only on L/£. 

A careful investigation of the model discussed by Chen and Dohm reveals that the power law contributions to the 
finite size corrections result from the interplay of two features of that model. The first is a sharp cutoff of fluctuations 
in momentum space and the second is the removal of the "remainder" term in ( |l.2j ) , which has the effect of introducing 
an interaction that cannot be periodic in reciprocal space. The combined effect of these two features is an effective 
interaction that falls off as a power law in the separation between degrees of freedom. This power-law interaction leads 
immediately to power-law contributions to the finite size corrections. The last result provides a natural explanation for 
the discrepancies between the finite-size effects that are predicted by infinite cut-off field-theoretical schemes [lj] 
and the finite size effects that arise in theories with a sharp cut-off ]TT| [l5^ . 

That a sharp cutoff is essential to the appearance of power-law corrections to the infinite system limit was noted 
by Chen and Dohm in |L2j , where an example of smooth cutoff effects is also presented. The both cases of sharp and 
smooth cut-off are clearly distinguished in |nj for systems with dimensionality d > 4. For d > 2 they also realized 
a close relation between a non-exponential large-distance behavior of the bulk correlation function (due to the sharp 
cut-off) and the power-law finite size behavior above T c fll5|| . Nevertheless, they do not consider in any of their articles 
the observed power-law finite-size corrections essentially as a result of an effective long-range interaction. 

In the present article we will deal only with the case when the hyperscaling holds, i.e. we will suppose 2 < d < 4, 
but definitely similar effects will be observed also above the upper critical dimension d u = 4, which case is a subject 
of an intensive discussion in the literature in the last time (see, e.g. jl9],^0) and the literature cited therein). 

It is certainly true that the degrees of freedom of a system of spins on a lattice is represented in reciprocal space 
by wave vectors confined to a Brillouin zone. However, the interactions between the spins have a Fourier-transformed 
form that is periodic in the zone. That is, the momentum-space version of the interaction is not purely quadratic, 
but rather reproduces itself as the wave number is shifted by a reciprocal lattice vector. Such interactions do not 
give rise to the non-scaling terms obtained by Chen and Dohm. An alternative source of the fluctuation cutoff is a 
Fermi surface. However, while the Fermi surface is a natural construct in the case of non-interacting Fermions, it is 
hard to imagine a property of any actual system that allows Bose excitations to be freely occupied in one regime of 
wave- vector space and that forbids the occupation of those excitations in an immediately adjoining neighborhood. No 
specific scenario leading to this behavior has yet been proposed, at least to our knowledge. 

To recapitulate, the two elements that are required to obtain the results of Chen and Dohm are, first, an interaction 
appropriate to a spin system in a continuum, and, second, a sharp cutoff in the fluctuation spectrum that mimics 
either the effect of a Fermi surface at T = 0, or a Brillouin zone. In Appendices [b|, |c] and [j5|, we explore the interplay 
of these two elements. We demonstrate that neither one alone suffices to to give rise to the effective long-range 
interaction that leads to a violation of finite size scaling. Appendix [a| supplies details of the analysis of Chen and 
Dohm that are central to the derivation of their predictions in terms of violation of finite size scaling. These details 
which are missing from their papers, are intended as an aid to those interested in a critical study of the basis for their 
results. The authors of this work wish to note that they are in full agreement with the mathematical conclusions 
drawn by Chen and Dohm. 

Regardless of the relevance of their particular model to either physical reality or the fictitious, but nevertheless phys- 
ically meaningful, world of simulations, Chen and Dohm have raised an interesting point. What can one reasonably 
expect to occur when interactions are intrinsically long-range, even if the range of the interactions is not sufficiently 
great to alter the asymptotic critical behavior of the system in question? Will these "subleading" long-range inter- 
actions inevitably alter the sorts of finite size effects that are naively expected to be present in finite systems with 
short-range interactions and periodic, or, indeed, arbitrary, boundary conditions? 

In this paper, we address this question in the context of an interaction whose Fourier transform allows for the kind 
of small-g expansion shown in Eq. (1.2). The term R(q), asymptotically smaller than the first two contributions to 
the right hand side of that equation when q is small, contains a component going as q a , where a/2 is noninteger, and 
a > 2. This means that the interaction is long range, but that 
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V(r)r p d d r 



(1.3) 



is finite for p < a. Alternatively, we imagine a V(f) going as r~ d ~ a for large values of r. Such an interaction is far 
from unphysical. In fact, the van der Waals interaction, which decays in three dimensions as r~ 6 , is consistent with 
(7 = 3. The explicit calculations presented in the article are for the case 2 < d < 4, 2<er<4 where d + a < 6. We 
find that an interaction of this form does, indeed, give rise to interesting modifications of the critical point behavior 
of a spin-like system. Those modifications include, but are not limited to, power-law finite system contributions that 
dominate the exponentially damped terms arising from a standard analysis of short range systems. In other words, 
in contrast to the infinite systems, where the subleading terms in the interaction are producing only corrections to 
the thermodynamic critical behavior, for finite systems they lead to dominant finite-size contributions in a given 
regime. The case a = d = 3 is special, in that there are logarithmic corrections to the nominally power law finite size 
corrections. 

The current understanding of the consequence of power law interactions with regard to critical point properties of 
an 0(n)-model spin system is fairly well-developed, but is not yet complete. Assuming that the interaction is of the 
form r~ d ~ a , then there are three regimes, depending on the magnitude of positive exponent a. First, if a < d/2, then 
the critical behavior is mean-field-like. On the other hand, if a > 2, the behavior of the system in the immediate 
vicinity of the critical point is as if the interactions are short-ranged. Clearly, when d > 4, the two regions overlap, 
and the critical behavior of the system is always expected to be as predicted by mean field theory. In fewer than 
four dimensions, the regime d/2 < a < 2 represents an intermediate case, in which the critical behavior is not 
necessarily dominated by either the mean field or the renormalized short-range interaction limits. This regime has 
been investigated for the 0{n) model in the context of a renormalization group-based expansion in 2a — d by Fisher, 
Ma and Nickel In the case of a one-dimensional Ising system, it has been argued that when a = 1, the critical 
properties are intimately related to those of the two-dimensional XY model pl[ , in that the appropriate version 
of the Renormalization Group is in the same generic class as the equations shown to apply to the latter model by 
Kosterlitz [^2j . The existence of a phase transition in this borderline case has been rigorously proven by Frolich and 
Spencer [g] . Later Imbrie and Newman ]23[ ] presented a rigorous proof of the existence in the system of a phase in 
which the two-point correlation function exhibits power-law decay with an exponent that varies continuously in a 
finite temperature range below the critical temperature. 

At this point in time, there seems to be no serious controversy on the analytical front. However, in a recently 
published paper Bayong and Diep [Q have utilized Monte Carlo simulations to investigate the critical properties of 
an continuous Ising system (i.e. the spin variables of the model can take any value between —1 and 1) in one and two 
dimensions with an interaction of the form l/r d+a . They seek to determine whether the behavior of this system is 
consistent with predictions based on Renormalization Group methods and other analytical approaches. While they 
are able to identify trends for the critical exponents that qualitatively follow the assertions of previous investigators, 
the values of the exponents are not in concert with those obtained in analytical investigations. This calls into question 
either basic assumptions with respect to the influence of long-range interactions on critical point properties, or on 
the validity of Monte Carlo methods, at least as exploited by the authors (e.g. for the case d = 1, a = 1.1 rigorous 
mathematical proof can be presented p5| for the absence of phase transition at finite temperature, while the authors 
determine the critical exponents around such a transition). Alternatively, there is the possibility that more care and 
experience is needed in the way in which long range interactions are discussed analytically. For example, will the 
long-range interactions studied by Bayong and Diep induce complications in the finite size corrections that require a 
more sophisticated approach than was undertaken by them? 

Additionally worthy of mention are the results of Luijten |2(| that indicate the existence of discrepancies between 
the Renormalization Group predictions for the behavior of the Binder cumulant B (for a definition of B see |27|) in 
a fully finite system with periodic boundary conditions at its bulk critical temperature and the numerical results for 
this quantity obtained by Monte Carlo simulations. Here, the focus is a system with leading power-law interaction. 
As a function of e = 2a — d Renormalization Group predicts (2^j2^] that B depends in a leading order on yfe, i.e. 
in the same way as in the case of short-range systems |T^ |, while numerical results [p6| suggest a linear dependence 
on e. The numerical simulations have been performed for a discrete Ising model with d — 1 and d — 2 and periodic 
boundary conditions. It is difficult to comment on the origin of this discrepancy — one is tempted to suggest a more 
careful analysis of the numerical data, despite the fact that the cluster algorithm |2j| used by Luijten is able to take 
into account the interaction of any spin with all the other spins in the system, including the infinite sequence of images 
under the periodic boundary conditions. In other words, no truncation of the interaction has been enforced. 

The outline of this paper is as follows. In Section ||, we outline the general features of the model to be studied: 
a mean spherical model in confined to a d-dimensional hypercube of length L per side, subject to periodic boundary 
conditions. The interaction contains a component with a diverging higher moment. The signature of this component 
is a term in the Fourier transform of the interaction going as q a with 2 < a < 4. As there is also a term that goes 
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as q 2 , the interaction is short-ranged, in that critical exponents are those associated with interactions having a finite 
second moment. Section [II contains a detailed analysis of the equation of state of this model, special attention being 
paid to the influence of the long-range portion of the interaction. It is found that an expansion in the strength of this 
contribution to the interactions between degrees of freedom in the system suffices to establish the key characteristics 
of the system — in particular, the finite size corrections to asymptotic critical behavior. This section provides all the 
mathematical detail needed to extract both the asymptotic critical behavior and the leading order corrections arising 
from non- leading-order long range interactions. In Section IV the results of Section III are utilized to discuss the 
dependence of the isothermal susceptibility of this system on both reduced temperature, t and system size, L in the 
important regimes surrounding the critical point. The case d = a = 3 is given special treatment, as in this special 
instance, which is relevant to a system in which van der Waals interactions play a role in ordering. Here, exponent 
"degeneracy" gives rise to logarithmic corrections to pure power law behavior. 

We find that one can easily envision situations in which corrections to scaling, in the form of contributions to the 
thermodynamics of a finite system that scale as L to a non-leading power are of a magnitude comparable to the 
putative leading order terms. This indicates that there are circumstances in which the analysis of simulation data 
must be undertaken with care. 



II. THE MODEL 



In this paper, we restrict our attention to the fully finite mean spherical model. We assume a d-dimensional 
hypercube of length L per side with periodic boundary conditions. The boundary conditions are consistent with the 
way in which one sets up a model for Monte Carlo investigations, in that the system is finite, but lacks physical 
boundaries. The degrees of freedom consist of a set of N localized spins with gaussian weight, and the Hamiltonian 
is given by 

H=- y £,v(R i -R )s i s j -hY,s i (2.1) 

where Ri is the position vector of the i th spin. The Fourier transform of the interaction, V (R4 — Rj^j > is assumed 
to possess the following low-q expansion 

v(q) = v Q + v 2 q 2 + v a q" + v A q 4 . (2.2) 



where 2 < a < 4, vq, v a > 0, and 1)2, V4 < 0. The term v a q a in (2.2) is associated with a contribution to the real-space 



version of the interaction going as \R\ — Rj\~ , as long as a is not an even integer (in the opposite case one will 
have in addition some logarithmic corrections). Note that the signs of the coefficients in the small q expansion are 
chosen so as they normally appear for interactions that decay in power laws with the distance between the interacting 
objects - molecules or spins (this easily can be checked, say, for the example of a one dimensional system with such 
an interaction; vq, V2, v a and V4 are tr-dependent — for simplicity of the notations this dependence is omitted here) . 
Furthermore, we suppose that v(q) — vq < if q 7^ 0, which reflects the fact that there are no competing interactions 
in the system we consider and that the only ground state of this system is the ferromagnetic one. Of course, it would 
be interesting to consider such systems — say a combination between antiferromagnetic short range and ferromagnet 
subleading long-range interactions, but this is out of the scope of the current article. 
The partition function of this system is equal to the multiple integral 



dsi--- j dsN exp [— f3H] (2.3) 
supplemented by the mean spherical condition 



N 



which can be enforced with the use of a "Lagrange multiplier" term going as A Xh=i s f the effective Hamiltonian, 
and thence into the partition function. The spherical model equation of state then takes the form 



E 



knT =N (2.5) 



A - (w + v 2 q 2 + v a q a + w 4 g 4 ) 
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The phase transition in this model occurs when the combin atio n (vq — A) /v2 = r takes on a value asymptotically close 
to zero. The difference between the equation of state in (2.5) and the standard mean spherical model con ditio n in 
short range systems lies in the addition of the term going as q a in the denominator on the left hand side of (2.5). In 
general, this term is taken to be negligible, but we will soon see that it leads to interesting effects. 

Because we are looking at a finite model, the sum over q in (2.5) is subject to restrictions. In particular, under 
the assumption of a periodically continued hypercubic system of length L per side, allowed values of q are of the 
form q — 2nn/ L, where ft is a vector with integer comp one nts. A number of methods have been developed for the 
evaluation of the kind of sum in the equation of state (^^). When there is no term going as q" , the sum is quite 
standard, and has been performed (neglecting the term proportional to fc 4 ) with the use of a number of techniques, 
including the Poisson sum formula and variants on the Ewald summation trick f30|| . The addition of the non-integral 
power of q into the denominator in (|2.5|) c omplicates matters a bit, but adaptations of the above methods to the 
case < a < 2 have proven effective [ |31| , |32| . Appendices [j| and |] outline such adaptations that incorporate also the 
case 2 < a < 4. They make use of contour integration tricks to "map" the summation onto the more conventional 
short-range one. In this paper, we make use of analytical methods based on the relationship between expressions 
central to the statistical mechanics of this system and well-known mathematical functions. 



III. THE EQUATION OF STATE 



The equation of state for the mean spherical model (for a comprehensive review on the results available for this 
is 



model see |10 



K 



Kr- 



+ WS(r,b,c,L\A), 



(3.1) 



where q = |q|, K = \v2\/ (ksT), b = v a /i)2 < 0, c = V4/V2 > 0, h is a properly normalized external magnetic field, and 

WS(r,b ) c,L\A) = j 1 J2 ' 



L d ^ r + q 2 + bq a + cq A ' 



(3.2) 



Here q is a vector with components qj = 2Trmj/L, mj — ±1, ±2, ■ • •, j = 1, 2, • • • , d, in the range —A < qj < A. The 
critical point of this system is given by K c (d, a, 6, c, A) = W%(0,b,c\A), where A = n/a, a being the lattice spacing 
for a lattice system (or A is the finite cutoff of the corresponding field theory), and 



W2(r,b,c\A) 



(27r) d 



d d 



1 



bq° 



cq^ 



(3.3) 



Wherever possible we will omit the contributions that are due to the term proportional to q 4 . Because of that we will 
omit c in th e re maind er o f the text in the symbols W% (r, b, c, L\A) and W?(r, b, c|A). The term proportional to q 4 is 
included in ( |3.2| ) and (3.3) in order only to assure that no artificial poles will exist in the denominator of W£(r, 6, L\A) 
and in the integrand of W%(r, 6j A) at large values of q (we recall that b < but c > and that the propositions we 
made for v(q) guarantee that there is no real root of the equation 1 + bq a ~ 2 + cq 2 = 0). 

We are interested in the b eha vior of the finite system close to or below the critic al t emper atur e K c , when r in the 
right-hand side of equation (3.1) is small (i.e. when < r < 1). As is clear from (|3.l|) and (3.3) the singularities in 
the behavior of W?(r,b\A) as a function of r, which in turn can be transformed as singularities with respect to the 
temperature, arise for small values of q. In what follows we will retain only those contributions to the behavior of 
the quantities involved that are associated with the effects of long-range fluctuations (i.e. q -C 1). Proceeding in this 
way, we obtain 



W° d (r, b, L\ A) ~ -L + Sl (d, r, 2| A) - £ 



rL d 



L d ^ ( r + q 2 ) 2 

q v 7 



+ S L (d, r, 2|A) - b(l + r^)S L (d, r, <r\A), 



where 



S L (d,r,a\A) = ^Yl 



q#0 



r + q 2 



(3.4) 



(3.5) 
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To analyze the finite-size behavior of Sl(<1, r, <r|A) we make use of the identity 



- — ■ — = / exp(-q 2 t)— iFt(a;a- p;-rt)dt, 

(r + q 2 ) a J r(a-p) 



(3.6) 



where a > p and \F\ is the confluent hypergeometric function. For a = 1 and p < 1 this identity further simplifies to 



2p poo 

-X"a= / exp[-( 9 a +r)t]rV(-P,-rt)*, 
r + q z Jo 

where t*(o, ar) is a single-valued analytic function of a and x, possessing no finite singularities 

x n 1 (-a;)" 

T(a + n+ 1) ~~ f(o) ^ (a + n)n! ' 



7*(a,«) = e 



X < oo. 



(3.7) 



(3.8) 



Both of the identities ( |3.6[ ) and (|3.7 i ) c an be proven by integrating by parts the corresponding series representations 
of iFi and 7*. In p3] the identity ( |3.6| ) has been used to analyze the finite-size behavior of 0(n) model with both a 
leading long-range interaction of the type \/ r d + 2 ~ 2a i Q _ > and a short-range interaction present in the system. 
With the help of this identity one obtains 



S L (d,r,2(p + l)\A) 



e- rt t- p i*(-p,-rt) 



q#0 



dt 



= IU(X, d\A) + L^-^I^rL 2 , d) + 

^ \ d+2p-2 

2^ 



M- 2 I P cut oS (r/A% 



where M = LA/(2w), and 



'bulk 



(r,d\A) 



e -rt t - Pl *(_ P) _ rt j 



d(2nAVi) 



/Ant 



dt, 



(3.9) 



(3.10) 



^scaling ( X ' ^3 



-*V-V(-p,-s*) 



E 

,k— — cc 



(4tt*) 



-d/2 



-1} dt, 



(3.11) 



and 



Jo 



e^ 47r2(1+;E)t erf(27rV / i)t 1 " P 7*(-p,-rf)^. 



(3.12) 



Taking into account the fact that j*(p,x) — > 1 when p — > (see Eq.(3.8)) it is clear that when cr = 2 all these 
expressions give the corresponding well-known results (see, e.g. fl2]| ) with only short-range term in the interaction. 
We will treat the bulk term separately. To derive the behavior of this term, it is not necessary to use the representation 
given here. In fact, it is relatively straightforward to derive the leading r dependence of the bulk term due to the 
existence of a subleading long-range term in the interaction. 

Furthermore, it is obvious that the term that contains a contribution due to the finite cut-off is not, in fact, correct, 
because the expansion we have utilized guarantees only that effects due to small q behavior are properly taken into 
account. In what follows we will simply ignore the precise form of this term. The "finite-size scaling term", insofar 
as it is due to long-wavelength contribution, is calculated exactly. We are able to conclude that the equation of state 
can be written in the form 



vL 



j2-djQ 

scaling 



(rL\d) 



d 



bL 4 - d - a (l + r ^-)I ( s ^ 2 {rL 2 ,d) + A-dcpcndcnt term. 



(3.13) 



(3.14) 



G 



In Appendix [g| we show that the A-dependent terms can be neglected for the purposes of the analysis carried out 
here. For the bulk term W% (r, &| A) for small r it can be shown by using the standard techniques (see Appendix ||) 
that 

wsirm = mom + ^r^sin^) ^ 2 " 1 

. j, n{d + a - 2) d/2-i+(o--2)/2 

2(4^) d / 2 r(d/2)sin(^(d + a)/2) 

+0(r (d+2(a-2))/2- 1)r)j (315) 

where 2<cr<4, 2<d<4 and d + cr < 6. 

Since d + a = 6 includes the especially important case d = cr = 3 we also present the corresponding result for that 
case 

mtr.MA) - >wia) + (4T) , / , r( ,; 2)sm(jrf/2) - j/2 - 



Inserting these expansions into the equation of state (3.14) we obtain 

2 

Xl = X sr (x) + bL 2 - a X lr (x) + h| + A-dependent terms + O (xL-^- d \x d / 2+a - 3 L- 2 ^- 2 A , (3.17) 



where x = rL 2 , Xl = [K -K C )L X I V ', x 2 = hL A / u /VK with f = l/(d-2), A = (d + 2)/(2(d - 2)) (see JToJ ) and 



A' r (x) 



Tr(rf + cr- 2) 



2(47r) d / 2 r(d/2) sin(vr(d + ct)/2) 



d/2-l + (cr-2)/2 



(3.19) 



To analyze the equation of state and the behavior of quantities such as the (reduced) magnetisation m — h/r, the 
susceptibility \ = 1/r JTo] the information required, in addition to that given above, is with respect to the asymptotics 



01 Scaling (M) &nd I 



ct-2)/2 
scaling 



(x, d) in different regions of the thermodynamic parameters. We will be interested in the 



behavior of m and x slightly above, in the region of, and below the critical temperature. 
The asymptotics of ^g Ca ii ng (^, d) are well known 



and 



where 



^scaling i. X i ^) 



1 

X 



dy/2 



^scaling ( x i ^) 



-^scaling (0; 



(3.20) 



(3.21) 



^scaling (0' ^) " 



dt 



dt 



k^O 

(47rt) 



-4?r 2 fc 2 t 



- (^t)- d/2 



k^O 



D 



(3.22) 



The corresponding asymptotics of Scaling v 2 ^ °0 are ( see Appendix |) 
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• Z 2aMng( !C ' d )- C * a! 2 ' Z^ 00 ' 



(3.23) 



where 



Co 



(i+p)4 1+ pr(i+ P + d/2) 



r d/2 



y 1 

Ld+2(p+l) 



r(-p) ^ fcd+2( P +i) 

fc#0 



and 



Scaling OM) - Valine (M) = 



1 



scaling 



r(i-p)y 
r(i-p)/ 



f/f 



(3.24) 



(3.25) 



Obviously the right hand sides are well defined for p < 1 and 2 < rf < 4 both around the lower and the upper limit of 
integration. We will denote this constant by D p . 



IV. FINITE SIZE EFFECTS AND THE SUSCEPTIBILITY 



Given the equation of state, we are now in a position to explore the behavior of various thermodynamic properties 
of the system with subdeading long range interactions. Here, we look at the susceptibility of such a system. We first 
consider the case 2 < d < 4 and 2 < a < 4. Furthermore, we assume d + a < 6. The scaling form of the equation of 
state is 



2 

51 = X sr (x) + bL 2 - a X lr (x) + ^| + A-dependent terms + o( xL^-^ , x d/2+<r-3 L -2(a-2) 

x 1 V 



Here, x = rL 2 , and the susceptibility \ is given by 



X = - 

r 



(4.1) 



(4.2) 



assuming that h = 0. We now analyze the behavior of \ for temperature, T, close the critical temperature, above and 
below, and for T in the immediate vicinity of the critical temperature, in that finite size rounding is evident. 

To that end, we need the asymptotics of X sr (x) and X lr (x) in the limits of large and small x. Making use of the 
asymptotics of / s ° caling (:r, d) and (1 + xd / ' dx)Il c&l{nE {x, d), we have 



U + A), 



30 

x -> 



and 



where 



and 



^ ' — ' D - - - - — >■ n 



(o— 2)/2j 



x — > 



(47r) d/2 r(d/2) sin (7rd/2) 



< 



Bd,a — 



7r(d + cr-2) 
2 (47r) rf/2 T(d/2) sin (7r(d + 2)/2) 

We begin with the case T = T c . The Eq. (fill) becomes 

= W 5, '(.t) + bL 2 ~ a X lr (x) + • • • . 



> 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 
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Let Xq be the solution of the equatio n X st (xq) = 0. Obviously, xq which is O(l), is a positive constant. Taking into 
account that a > 2 and solving Eq. (4/7) iteratively, we obtain 

I ~ 1 + bL 2 -°X lr (x )/ {x 2 X' sr (x )) (4.8) 

X Xo 

where X' sr (xo) is the derivative of X sr (x) at x — xq. Recalling that \ = V r an d x = r ^ 2 ; one immediately obtains 
from (gj) 

X ~ x^L 2 + bL 2 -^X lr (x )/ (x 2 X' sr (x )) (4.9) 

It is clear that in a Monte Carlo simulation if neither a nor L is particularly large, then the correction terms in ( |4.S| ). 
which go as i 2 ^" 7-2 ) might well be as large, numerically, as the leading order terms which scale as L 2 , depending, of 
course, on the values of b, X lr (xo) and X' sr (xo)- 

Let us now consider the case in which T is fixed close to, but also above T c and L — > oo. Then, x± — > — oo, and, 
t akin g into account the corresponding asymptotic behavior of X sr and X lr for x — > oo (rL 2 3> 1) we can rewrite Eq. 
(4.1) in the following form 

K - K c ~ Ay' 2 ' 1 + bB d ^r d ' 2 - 1+a / 2 - 1 + bC p L- (d +°\- 2 (4.10) 

Solving this equation iteratively, we obtain 

>-2)/2 



X - Xo < 1 + bjXo 



B d,a . n ( 2 \( d +°-)/ 2 
V L-(<t-2)/2 \X0 L ) 



A d 

where xo is the susceptibility of the corresponding infinite system with short-range interactions only, i.e. 



(4.11) 



with 



The above solution is valid when rL 2 ^S> 1, i.e. L 2 ^> xo- Note that the dominant finite-size corrections to the 
behavior of the total susceptibility are of order L~( d+a \ That is, they are not exponentially small, nor are they 
cutoff-dependent. The existence of corrections of this type in the case of leading-order long-range interactions is well 
known. First they have been derived in the framework of the spherical model [ p4| , ^5| . Analogous is also the behavior 
of the 0(n) model within e = 2a — d expansion j2{| (at least up to the first order in e). 

Finally, let us consider the case T < T c . Then, X\ — > oo, which leads to x — > 0. Eq. ( |4.l|) becomes 

K - K c ~ -ij + D Q L 2 - d + bL*-^D (a _ 2)/2 + ^- 2 . (4.14) 

In the absence of an external field, the iterative solution of the above equation yields 

X ^(K- K c ) L d - D L 2 - bD (a _ 2)/2 L 4 -° (4.15) 

Now, we turn to the case d + (7 = 6(2<<i<4, 2 < a < 4). This is especially a propos, in light of the fact that the 
van der Waals interaction in three dimensions leads to a contribution in which d = a = 3. In this case, it is necessary 
to take into account the special form of X lr : 

d \ T (4-d)/2, 



X lr {x) = 2Bx\nL - Bxlnx + \1 + x— j (x, d), (4.16) 
where 

B= j£ . (4.17) 

(47r) d/2 T{d/2) 

The first term is responsible for the leading-order finite-size corrections that are due to the subleading long-range part 
of the interaction. 

Proceeding as in the case d + a < 6, it is readily demonstrated that 







a) For T = T c : 



(4.18) 



The correction term is obviously important in the analysis of Monte Carlo data, 
b) For T > T c : 

B 1 



X-Xo 



(4.19) 



c) For T < T c : 



In this region, the corrections to bulk behavior due to long-range corrections play no role, and the solution 
remains unaltered. 

V. CONCLUSIONS 

In this paper, we have reported the results of an investigation into the critical point properties of a finite spherical 
model in which interactions contain a component that is long-range, but insufficiently so to alter the asymptotic 
singularities of its thermodynamics — in particular the critical exponents. One can envision interactions decaying as 
r —(d+a) ^ w j lere d ig dimensionality of the system and a > 2. An important example is van der Waals interaction, 
which decays in three dimensions as r~ 6 that is consistent with a = 3. The finite system that we consider is subject 
to periodic boundary conditions, and, thus, provides a model for the sorts of systems that are studied in computer 



simulations. This investigation was stimulated by recent work of Chen and Dohm, |1T 13j in which a combination of 
a spin-spin interaction truncated in momentum-space and a sharp momentum-space cutoff on fluctuations gives rise 
to an effectively long-range interaction. 



In the critical region we find that the susceptibility of the finite system \{t, h; L) is of the form (see Eq. (4.1), ([1.2]), 



dyD) 



or, equivalently, 



X (t,h;L) = L^ v Y(x 1 ,x 2 ,bL 2 -"), (5.1) 



X (t, h; L) = Lr/ V [Y sr {x u x 2 ) + bL 2 -°Y lr {x u x 2 )\ , (5.2) 



where x\ = citL 1 ^, x 2 = c 2 hL A /", and Y, Y sr and Y lr are universal functions. The quantities c\, c 2 and b are 
nonuniversal constants. Note, that the above structure of the finite-size scaling function in systems with sublcading 
long-range interactions is different from the corresponding one for systems with essentially finite range of interaction 
]^To|. The new length scale which is involved does not lead to corrections of the finite-size scaling picture known 
before, but leads, see below, to leading finite-size contributions above T c . 

In the range of parameters, for which tL 1 l v — 0(1), and also a bit below the critical point, where tL 1 / u — > — oo, 
the long-range contributions represented by Y lr are merely corrections to the leading finite-size behavior. Somewhat 
more interestingly — and of greater practical significance — there are also corrections to the dependence on system size 
of singular thermodynamic properties at the critical point that can conceivably cloud the numerical analysis of Monte 
Carlo data, in that the corrections, while less important in an asymptotically large system, may be of the same order 
of magnitude in systems that are a realizable size (see Eqs. (4.9), (4.15) and ( 4.18| )). The studies reported here ought 



to provide, at the very least, a conceptual basis for the critical evaluation of Monte Carlo results. 

On the other hand, in the high-temperature, unordered phase, where tL 1 ^ — > oo, we find that the long-range portion 
of the interaction between spin degrees of freedom gives rise to contributions of the order of bL~^ d+a ' that swamp 
the exponentially small terms that are expected to characterize the signature of finite size in systems with periodic 
boundary conditions and short range interactions. In other words the subleading long-range part of the interaction 
gives rise to a dominant finite-size dependence in this regime. This is entirely consistent with the inherent long-range 
correlations that attend long-range interactions, but it goes beyond the standard finite-size scaling formulation. More 
explicitly, one obtains Y sr (xi : Q) ~ exp(— const, x^), while 
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(5.3) 



when x\ — > oo. This asymptotic follows from the requirement the finite-size corrections to be of the order of in 
this regime, which is to be expected on general grounds and is supported from the exis ting both exact and perturbative 
results for models with leading long-range interaction included [ p4|j35|j2^] . Note that (5.3) implies for the temperature 
dependence of this corrections that 



X (t, h; L) - X (t, h; oo) ~ t -{d+2)u--y L -{d+a) t j\/u 



(5.4) 



This prediction is in full agreement with Eq. (4.11). Obviously, the existence of such power-law finite-size dependent 
dominant terms above T c is of crucial significance in the analysis of the Monte Carlo data from simulation of such 
systems. 

Finally, it is worth noting that the system considered here is equivalent to an 0(n) dimensional vector spin model in 
the limit n — ► oo [(To). Because of the spin- wave excitations, the bulk correlation length of such a system is identically 
infinite below T c (for any 0(n), n ^ 1 model). As a result the correlations decay in a power-law in this regime. 
The direct spin-spin interaction decays faster there, and that is why for T < T c we obtain no essential finite-size 
contributions due to the subleading term in the Fourier transform of the interaction. The situation is different in 
Ising-like systems. There below T c the correlation length is finite, the correlations decay exponentially fast in a 
system with only a q 2 term in the Fourier transform of the interaction. Since, when the interaction is long-ranged the 
correlations cannot decay faster than the corresponding direct spin-spin interaction, one should expect modifications 
of scaling of the type presented in the current article for T > T c to be necessary for Ising- like systems also below T c . 

In addition, in the current article we have proven the approximation formula (see Eq. (A14)) 



M-l 

£ « 

n=-M 




erf (Vi^Mj +]Te- fe /4t 



-Mte 



—AirtM 



(5.5) 



which is of a bit more general mathematical interest and which also might be useful in a lot of studies of finite-size 
effects by exact or perturbative methods. 
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APPENDIX A: DERIVATION OF THE CENTRAL RESULT OF CHEN AND DOHM 

The expression in which the size dependence of the statistical mechanics appears in the papers of Chen and Dohm 
p] |l3[ is given by the difference between a sum over the set of allowed wave vectors in a hypercubic system having 
a linear extent L in every direction and the integral for such a system in the limit L = oo. It is assumed that both 
the sum and the integral are taken over a region of wave-vector space that is also a d-dimensional cube. The system 
under consideration is subject to periodic boundary conditions in all dimensions. This paper addresses the question 
of the source of the violation of scaling found by Chen and Dohm by focusing on the effective long-range nature of 
the interactions that are generated by the combination features assumed to hold for the system considered by them. 
However, for the reader interested in looking at their papers on the subject we provide here details of the derivation 
of the terms in the equation of state that fall off as a power in the size, L, of the system. 

The quantity from which the power-law finite-size corrections arise is the difference between a lattice sum over 
wave-vectors, k and the integral to which it is equal in the thermodynamic limit. This difference, which is introduced, 
for instance, in Eq. (5) of p3| , is given by 



Ax (A, X - 1 )=I(A,X- 1 ) 



Sl (A,x _1 ) , 



(Al) 



where 




(A2) 



and 



5 L (A, x - 1 )=L- d E 



1 



(A3) 



n^0 X 



+ 



2tt 
L 



with 



n,i e [-M, M), i = 1, • • • d, M = LA/(2tt). 



(A4) 
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We start by looking at the term 7 (A, x This term yields readily to analysis. Making use of the identity 



and 



1 

2^ 



f°° 1 
/ exp (— ax) = — , Re a > 0, 
Jo a 

erf [Ay/t\ 



exp [—ak ) dk 



-A 



the function I (A, \ 1 ) can be rewritten in the form 

f°° r r i ~i 1 

I (A, x" 1 ) = J (4nt)- d/2 exp [-ix" 1 ] {erf [AViJ } 

We now turn to the expression for the sum, Sl(A, x _1 ). This sum is given by 



X _1 L 2 + 47r 2 n 2 



= L -(d-2) 



M-l 



1 



\n=-M 



dt-L- d —-. 

X 



We now write 



JVf-l 



S(t,M)= E 



=-M 



and apply the Poisson summation formula 

b 



£/(")= E / e^ fc "/(n)dn+i[/(a) + /(6)] 



/c— — oo 



This yields 



00 „M 

S(i,M) = E / e 

k=-oo J - M 

1 ^ /.JM 

= -= erf ( VA^HM ) + 2 V / g- 4 ^ 2 ™ 2 * CO s 2tt 



-4:7rtn -\-2nikn 



kn dn 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



(All) 



If the only term retained in our result for S(t,M) were the first one on the last line of ( All ), then there would be 
perfect cancellation between the sum and the integral. The difference between the two results from the second term. 
We now construct an asymptotic expansion for that term. We have 



°" pM °" f PC 

2E / e- i7r2n2t cos2irkn dn = 2^2\ 
k=i J - M k=i 



; -47r n-t CQs27Tkn dn _ 2 



' cos 2irkn dn 



M 



k— — c 
oo 



k— — oo 
oo 



k=— oo 



k 2 /it 


- 29? 


\ 






Jm 






' poo 


k 2 /it 


- 23? 


Jo 








k 2 /it 


- 8Mte~ iw 



e- 4 * tn+27tikn dn 

e -iTT 2 t(M+n') 2 +2TTik(M+n') dn l 
1 



■LM 2 



k 2 + 16ir 2 t 2 M 2 



(A12) 
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We obtain the last line in ( A12| ) by retaining only those terms in the exponent in the next-to-last line that are linear 
in the integration variable n' . 

It is now fairly straightforward to demonstrate that the combination k 2 + l&n 2 t 2 M 2 in the last line of (A12) can 
be replaced by k 2 with no loss of accuracy in the evaluation of the sum. Making use of the result 



00 i 2 

El 7T 
12 ~~ ~~Q 



we end up with 



M-l 

£ - 

n=-M 



k=l 
1 

4tt 2 



(A13) 



erf (VAttHm) + ^ 



,-k 2 /it 



k^O 



-Mte 



-A-k tM 



(A14) 



The remainder of the calculation involves the insertion of the above results into the expression ( A8). The key term 
arises from a cross-term in the expansion of of the d th power of ( A14 ). That term is 



UirhM (Airt)- (d - 1)/2 [erf (V^ 



HM 



(d-l) 



iTT A tM A 



(A15) 



The remainder of the calculation involve scaling the system size, L out of the integral over t in ( AS ) . To recover the 
form exhibited in E2], on makes use of the equality 



(47rf) 



-(d-l)/2 



tM 



d-l 



1 

2~t7 



M 



-M 



-tk- 



dk 



d-l 



and of the Jacobi identity 



k^O 



k 2 /4t _ 



/4^^e- 4rft 



(A16) 



(A17) 



which leads to the end-result 

S L (A, x _1 ) = I (A, X" 1 ) - L-^h (L 2 x^) - K d - 2 ai (d, x~ 1 A -2 ) (ALy 2 + ■■■ 
which leads to 

A x (A, x- 1 ) = L-^h (L'x- 1 ) + A d - 2 ai (d, X -'A- 2 ) (AL)- 2 

with the functions I\ and a% as given by Eqs. (7) and (8) in (l3[ , respectively. Note that h(x) = /s Caling (x, d), and 
a 1 (d, 2 ;) = -/ c utoff (x)/(47r 2 ). 



(A18) 
(A19) 



APPENDIX B: THE COMBINED INFLUENCE OF AN INTERACTION GOING AS Q 2 AND A SHARP 

CUTOFF IN Q-SPACE 

Consider an interaction tha t Fo urier transforms to q 2 exactly. That is, imagine that the Fourier transform on the 



interaction is as given by Eq. (1.2) with R(q) — 0. The first term in this expansion gives rise to a real-space interaction 
that is entirely local. We thus focus on the term that goes as q 2 . For simplicity, we start by restricting our attention 
to a one-dimensional system. If the interacting spins reside on a lattice with unit lattice spacing the form of this 
interaction in real space is given by 



d 2 
dn 2 



q 2 e^dq = ~—J e tqn dq 

dn n 
^7T 3 n = 



3 cos™ , Q (Bl) 
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As n is an integer, the real space interaction decays as a modulated power law. Such an interaction has a different 
range than one in which the power law is "pure," in that there is no alternation in the sign of the interaction as the 
distance between the spins increases. However, this system displays long-range correlations, which manifest themselves 
in the size-dependence of thermodynamic quantities. 

One can assess the impact of this interaction on, say, the free energy by separating it into a short-range piece and 
a long range one. At one extreme, one can take the short range piece to be the n = part of the interaction. If 
one retained that term and discarded the long-range portion of the interaction, one would be left with a system in 
which there are no correlations between degrees of freedom at any non-zero separation. Suppose we start with this 
approximation, which is not too far from reality at very high temperatures. Then, we discover how the long-range 
portion of the interaction influences a system with periodic boundary conditions with the use of the perturbation 
expansion. Writing the Hamiltonian of this system in the form 



H = H + H h 



(B2) 



where Hq is the short-range portion of the interaction, cx sf , and H\ T is the long-range portion, we write the free 
energy as follows 



T 



-In 


E e ~ w 




.{Si} 


-In 


E e ~ H °~ 




{Si} 


-In 


E^° 




.{Si} 


F + 


(H h ) + ■■■ 



(B3) 



The first term on the last line of (B3) is the free energy of the system with short-range interactions only. The second 
term is the result of expanding the free energy to first order in the long-range portion of the interaction. This second 
term takes the form 



V(n)(s 



(B4) 



where, V(n) is the real space version of the long-range interaction, as given in Eq. (Bl). 

Let's imagine the case of a system with periodic boundary conditions. Such a system is equally well represented 
by an infinite set of duplicates of a the finite system. In this case, the correlation function (s m s m+n ) is equal to zero 
unless n — or n = ±kL, where k is an integer and L is the size of the system. We take L to be an even integer. One 
then obtains for the influence of the long range interaction on the free energy of the system 



oo 

,2 \ \ " 



1 



^ k 2 L 2 
fc=i 



(B5) 



Note that this influence goes as a power law in the size of the periodically continued one-dimensional system. It is 
interesting to note that the detailed dependence on the size of the system is different from the above when L, the 
system's size in terms of the distance between sites, is an odd integer. 

In three dimensions, the one-dimensional Brilouin zone is replaced by a cubic zone in reciprocal space, and the 
interaction in real space is given by 



dq x / dq y 



dqJql + ql + ql) e^ H 



(B6) 



It is straightforward to see that the interaction consists of three contributions, each long range in one direction and 
extremely short range in the two others. For example 



dq x / dq y / dq z q 2 x e q ' n = V{n x )5 ny 8 nz 



(B7) 

The overall interaction is thus long range, but highly anisotropic. That is, a given spin interacts with spins arbitrarily 
far away, but only with spins separated from it by a displacement vector that points entirely along the x, y or z axis. 
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APPENDIX C: THE EFFECT OF TRUNCATION OF THE FOURIER TRANSFORM OF THE 

INTERACTION 



In Appendix [B|, it is established that an interaction whose Fourier transform is truncated at the quadratic term, 
coupled with a sharp cutoff in Fourier space, has, in real space, a modulated power-law tail. In this and the following 
appendix, we demonstrate that both the truncation and a sharp cutoff are required for that long-range behavior to 
be manifested. Here, we investigate the effect of truncation only. As in the previous appendix, we focus on one 
dimension. There is every reason to believe that the our conclusions are unaltered in higher dimensionality. 

We will limit our discussion to the sum over wave vectors entering into the spherical model equation of state. First, 
we evaluate that sum for a lattice system with nearest neighbor interactions. In that case, the Fourier transform is of 
the form cos q, where q is the wave vector. We assume unit lattice spacing. The sum of interest has the form 



i=l 

where 



N 1 

E — - — ( C1 ) 

a — cos qi 



*i = 2 4 (C2) 



N 

This sum can be recast as a contour integral. Write 

cosq = Uz+^\ (C3) 

where z lies on the unit circle. The z's appearing the sum in ( |C1| ) are the N values of the N th root of unity. The 
contour integral version of the sum is 

1 Nz N ~ l 

w TT ~M dz (C4) 

a- \{z + l/z) z N -l y ' 

The contour is actually a set of TV contours, each going counter-clockwise about one of the N roots of z N - 1. To 
express the temperature dependence of the sum, we replace a by 1 + t. The roots of a — k (z + -) then lie on the 
real z-axis, inside and outside the unit circle. If we call them n and r^, where r\ is the root that lies outside of the 
unit circle, then 



n = 1 + 1 + V* 2 + 2t (C5) 

while ri = l/r\. The result of the deformation of the contour is a sum of two terms, one associated with a contour 
encircling each of the the two roots r\ and r-i . The end result is 

When t is small, 

ri w (C7) 

The corrections to the infinite system limit go as e~^^ N . 

For comparison, we now approximate cosine functions by the first two terms in the expansion in their arguments. 
The sum in (plf) is replaced by 



^ 1 1 f N 2 



^ t + & 2 2m J c 2tt 2 M^t + z 2 tan wz 

i=— oo N z L 2n' ! 



dz (C8) 



where the contour of integration is a set of contours, each circling counterclockwise about the zeros of tan7rz, which 
lie on the real axis at integer values of z. These contours can be deformed into two, one going around the zero of 
M-st + z 2 in the top half plane and the other going around the zero of that function of z in the bottom half plane. 
The evaluation of residues leaves us with the final result 
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N 



(C9) 



The expressions (|C6| ) and (C9) have the same limiting values at nonzero t, N — ► oo and at large N as i — > 0. In fact, 
one can rewrite (|Cq ) as 



JV(1 + V2i) 



1 



— 1 tanh Ny/t/2 y/2t + t tanhTV, 



(CIO) 



The exponentially small finite size corrections are, in leading order, identical. 

Thus, a truncation of the Fourier transform of the interaction potential does not influence the asymptotic form of 
the finite size corrections, if there is not also a cutoff in momentum space. Of course, some sort of cutoff is required 
in two or more dimensions in order that the sum with truncated interaction potential does not suffer an ultraviolet 
(large q) divergence 



APPENDIX D: THE INFLUENCE OF A SOFT CUTOFF ON THE RANGE OF THE INTERACTION 

The long-range interactions and correlations derived in Appendix [b| are the consequences of both a truncation in 
the g-space expansion of the interaction and the existence of a sharp cutoff in momentum space. To see how the 
"roundedness" of the cutoff changes the range of interactions, we modify the sum over q in the one dimensional space 
by introducing a cutoff function, having the form 

C(q,Q) = TT ^ TJW (Dl) 

In the limit Q — > 0, this cutoff approaches a step function. For finite Q, the gradual nature of the cutoff causes the 
interaction to be intrinsically short-ranged. The real-space version of the truncated interaction is, in one dimension, 

f°° (a + bq 2 )e iqx , 
V(x)= I l + JJ 1)/Q2 dq (D2) 



The behavior of the interaction can be extracted by distorting the contour of integration in (D2) so that it encloses 
the poles of the cutoff function, which are at the locations in the complex q plane at which 

l±*Q 2 fn + ^7rV 2 «l±i7rQ 2 ^i (D3) 

where the approximate equality holds if Q is small and the arbitrary integer n is not too large. The residue at such 
poles has the ^-dependence e~® ~" . This exponential decay of interactions, and hence of "intrinsic" correlations, 
will not give rise to effects interfering with the finite size corrections that go as e~ L ^. This can be seen by, first, 
repeating the analysis of Appendix O. Alternatively, one can look at the sum 



1 



^ r + g 2 1 + e (<?-i)/Q" 

9 



(D4) 



where the allowed values of q are 



27T71 ,„ , 

9 = -T (D5) 



Again, L is the size of the system in units of the lattice spacing between spins. The sum in (D4) is evaluated with 
the use of the Poisson sum formula: 

E/(«) = E / dqf(q)e lmqL dq (D6) 

q m— — oc 

The m = contribution to the sum in question is just the infinite system limit of the sum in the equation of state. 
Finite size corrections arise from the m ^ terms in that sum. Such terms have the form 
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imqL 



r + q 2 l + e (« 2 -i)/Q ; 



-dq 



(D7) 



We evaluate the integral in the same way as we did in the case of the interaction, except that here there is an 



additional pole in the complex g-plane, at the root of the denominator in l/(r 
a contribution going as e 



q 2 ). The residue at this pole yields 



as e 



ms/rh _ e mL/^ rpj^ j-ggj^gg f p \ es f the cutoff function give rise to terms going 
when Q is small and n is not too large. As the critical point is approached, r becomes small, the 



contribution going as e ^ L dominates all others, and finite size scaling in the expected form is recovered. 



APPENDIX E: ALTERNATIVE APPROACH TO THE ANALYSIS OF THE EQUATION OF STATE 



The principal influence of the subleading long-range contribution to the spin-spin interaction is obtained by ex- 
panding to first order in that interaction. In the case of the equation of state, the correction term, exhibited on the 
first line of Eq. (3.4), is 



E 



r + k 2 ) 2 



(El) 



The exponent a is greater than one. We start by making use of the contour integration identity 



(k 2 



1 

2ni 



7/2 



1 



(z + r y 



k 2 



dz 



(E2) 



where the integral is around the contour that encircles the pole of the integrand at z = k 2 . Figure [l] shows the contour 
over which the integration is performed. The integrand in ( |E2| ) has, in addition to the abovementioned pole, a branch 
point at z = 0, which results from the term z a ' 2 , assuming that the quantity a/2 is not an integer. There is also a 
double pole at z = —r. The integral in (E2) is evaluated by deforming the contour so that it surrounds the branch cut 
from z = to z = — oo, with a special accommodation at the pole at z = —r. Figure || is a picture of the deformed 
contour. The integral over this new contour has the form 



sin (T7r /2 



W 2 1 d 

7 \2 ,i,2 dz + cos (W 2 ) T 

(r — z) z z + k A dz 



W 2 



z + k 2 



(E3) 



In this expression, the integration variable z has been replaced by — z. The integration in the first term in (E3) is 
understood to be in the form of a principal parts integral. Such an integration combined with the second term has the 
effect of removing the non-integrable singularity at z = r. A new form for (E3) is obtained with the use of integration 
by parts. The end result is an integral of the form 



cos 7TCT/2 



dr 



r a/2 

+ k 2 



sin7r<7/2 



— z dz 



W2 



k 2 



dz 



(E4) 



The integral in (E4) is understood as a principal parts integral. 

The next step is to perform the sum over k. This sum is fairly straightforward, in that it is the one encountered in 
studies of finite systems with short range interactions p8| , p7{ . That has been done previously. In d dimensions, the 
result of the summation is given by 



Y — 

^ z + k 2 



L 2 
-in 



l 

v_ 

47T ./ | 



1] dt 



(2n)° 



d d k 
k 2 + z 



(E5) 



where 



Q(t)= e ~ 



n— — oo 



(E6) 
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If we extract the infinite system term — the last term on the right — from the right hand side of Eq. flE5| ) , we are left 
with a function that has the general form 



L 2 g(zL 2 ) 



(E7) 



In addition, this function decays exponentially with large values of zL 2 . In fact, it goes as e~ zL / 4lr . For small values 
of z, the function is dominated by the term 1/zL 2 . 

Figure |^ is a graph of the sum, with the infinite system term removed. Given the general form of this function, we 
can write the expression for the correction to the leading order contribution to the equation of state. It has the form 



cos(<77r/2)^ [r^ 2 L 2 g(rL 2 



sin cr7r/ 2 



z dz 



z a l 2 L 2 g(zL 2 



dz 



(E8) 



The expression above for the correction to the equation of state can be shown to have the expected properties in 
various regimes. For example, if rL 2 is large, then the leading order behavior arises from th e se cond term in (|E^). If 
we rescale the variable of integration by making the replacement z — > zj L 2 , the integral in (E8) becomes 



Sill (77T 



4-a 



"L 2 — z dz 



z a/2 g{z) 



dz 



-4-a 



G(rL 2 ) 



(E9) 



At large values of rL 2 , the integrand can be expanded in inverse powers of that combination. The lowest order term, 
going as 1/ rL 2 , can be shown to integrate to zero. The next order term is 



sin an L a 

2 

n r z 
Finally, use of the identity 



dz 



W 2 , 



dz = — 



sin an L 



7T 



W 2 



g(z) 



dz 



r p-i 



-dx = — n cot pn < p < 1 



(E10) 



(Ell) 



allows us to show that the expression in Eq. ( |E8| ) for the correction to the equation-of-state sum does not give rise 
to any terms going as a fractional power of r in the limit r — > 0. This is consistent with the expectations one has for 
the limiting behavior of a finite system, and with the analysis in Sections III and |v|. 



APPENDIX F: THE EQUATION OF STATE IN GENERAL 



While an analysis sufficient for our purposes can be carried out by focusing entirely on the first order effect of the 
long-range component of the interaction to, say, the equation of state, it is also possible to write down an expression 
for the entire equation of state for the system with a subleading long-range interaction. We start with the identity 



1 



ak 2 + bk G + ck 4 



1 

2ni 



1 



1 



+ az + bz a / 2 + cz 2 z — k' 



-dz 



(Fl) 



Where the closed integration contour is about the pole in the integrand at z — k 2 . See Figure [ij 

The next step is to distort the contour so that it wraps around the two poles of l/(r + az + bz a / 2 + cz 2 ), and 
around the branch cut on the negative z-axis. These poles exist for the range of a considered here: 4 > pa > 2. The 
distorted contour is as depicted in Figure ^. In the case of the two roots in the right half of the z-plane, there is a 
contribution from the residue proportional to the inverse of the derivative with respect to z of r + az + bz a I 2 + cz 2 . 
For the contour integral, because the direction of integration is different on the two sides of the branch cut. We take 



the difference between l/(r 



E 



b(ze 



■ CZ 



l/(r- 



bz a l 2 cos an j2 ± ibz a/2 sin an/2) 



ak 2 + bk a + ck 4 



1 



a + b{a/2)z a 



a/2-1 



2czo 



k 2 



bz a l 2 sin an j2 



(r - az + cz 2 + bz a l 2 cos an/2) + b 2 z a sin 2 an/2 k 2 + z 



-dz 



(F2) 
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In Eq. (F2), the quantity zq is the solution of the equation 



,<r/2 



= 1 



(F3) 



There are actually two solutions to (F3). When a Z, 2, they lie just above and below the negative real z-axis. In fact, 
as r — ► 0, zq — > —r/y/a± iS, with i«r. We will assume that the zq that enters into Eq. (F2) is the solution lying 
in the upper h alf of the complex z-plane. The sum over k in ( |F2] ) is carried out in the standard way, the result being 
given by Eq. (E5). 

While the resulting expression for the sum over k in the equation of state is correct to all orders in the amplitude 
of the singular contribution to the interaction, the expression presents significantly greater challenges to the theorist, 
and all important effects are recovered from the first order term in the expansion with respect to the long-range 
interaction. 



APPENDIX G: ESTIMATE OF THE LEADING A DEPENDENCE OF THE FINITE-SIZE TERM 



We are interested in determining the leading A dependence of the sum 

W%(r,b,L,c\A) = ± i J2 ' 



L d ^ r + q 2 + bq a + cq 4 ' 



(Gl) 



We will assume a sharp cutoff, i.e q is a vector with components qj = 2irmj/ L, rrij = ±1, ±2, • • •, j = 1, 2, • • • , d, in the 
range —A < qj < A. The leading A-dependence of the sum arises from terms with large q. One can, therefore, omit 
the q 2 and q a contributions. Let us denote the leading A-dependent term of W%(r, b, c, L\A) by AW%(r,b,c,L\A). 
Then 



M-l 



M-l 



aw: 



rf >Ac,z,|A) 1 E ■•• E YTT^ 2 



y ' y mi=-M m d =-M v 

where y = (r /b)(L/2ir) i and m = {rnj H h m 2 d ) 1/2 . Using the identity jlO) 

= / dt e- zt t a - Y E a a (-t a ), 

where E a>a {x) are the Mittag-Leffler functions, the above expression can be rewritten in the form 

I 



(G2) 



(G3) 



AW%(r,b,c,L\A) 



dxxE%fl(—% 



M-l 

E 



-xm /y/y 



=-M 



(2irybyL*- 2 J 

Taking into account that E2^{z) = sinhz/y^ |l(| with the help of Eqs. ( |A14| ) and ( A15| ) we obtain 



AIF d CT (r,6,c,L|A) 



487r( 9 - d )/ 2 



ML 



d-4 



dxx (5-d)/2 Mxy/y) d ~ x „- . > / - < 



x^/y 



[erf(MVa)]' 



(G4) 



(G5) 



Since M 3> 1 in the above equation, the only significant contributions are those stemming from small x. After taking 
into account that M 2 ^> y 2 ^ a , and, so, one has \im x _, sm(x^/y)/(x^/y) = 1, we are led to 



AW%(r,b,c,L\A) = 0{L- {d -^M d ^) = 0(M~ 2 A d ~ 4 ). 



(G6) 



The contribution in the equation of state that does not depend on A, but is due to the long-range character of the 
interaction is of order 0(L 4_d ~ CT ). Elementary checks reveal that when d < 4, d + a < 6, L^>1 and A ^> 1 one has 
j~6—d—(r > A <i-6 5 w hence, 



AW%(r, b, c,L|A) < L 



4-d-CT 



(G7) 



Thus, up to the order at which the results are presented in this article, those results will not be influenced by the 
(nonuniversal) A-dependent corrections. 
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APPENDIX H: DERIVATION OF THE LEADING ASYMPTOTICS OF THE BULK TERM W£{R,B\A) 

As we are interested only in the contribution stemming from small q, the expansion below is justified and one 
obtains 



1 



'-A r + q 2 + bq IJ + cq 4 
= W^MK)-r^ d jy q ...J_ 



A q 2 (r + q 2 ) 



„CT-2 



+br——-T / d q ■ ■ ■ , 

(2nyJ_ A J- A (r + q 2 ) 2 



i^) d J-A H J-Ar + q 2 
i rA ,A 2(<r-2) 

-^7^/ A"--/ 4t"5«+-- (HI) 



The nonanalitycity in the behavior of all this function arises entirely from small-g contributions in the integrals. We 
sphcricalize the region of integration and extend the limits of integration from zero to infinity. If the corresponding 
integral diverges after such a procedure, we first differentiate the requisite number of times with respect to r, perform 
a replacement of the limits of integration in the first derivative that does not diverge, calculate the leading r behavior, 
and, finally, integrate the required number of times with respect to r. Performing this procedure we obtain (d + a < 6, 
2<rf<4, 2<cr<4) 



WS(r,b,c\A) = WS(0,b,c\A) 
+ 



d/2-l 



(4ir) d / 2 r(d/2) sin(ird/2) 



ir(d + (J-A) „d/2-l+(a-2)/2 



+fe 2(47r)' i / 2 r(d/2)sin(^(d + ( T)/2) r 



d/2-l+(<r-2)/2 



+b (An) d / 2 T(d/2) sin(Tr(d + a) /2) T 

+ 0(r (d+2( CT -2))/2-l jr)j (H2) 



wherefrom we are able to obtain the result given in the main body of the article. 
In a similar way, one can treat the case d + a = 6 with 2 < d < 4, 2<er<4. 



APPENDIX I: DERIVATIONS OF THE ASYMPTOTIC FORM OF THE NON-LEADING LONG-RANGE 

CORRECTION TERM 

In this appendix, we provide details of the calculation leading to the asymtotic form of the expression 
(1 + xd/dx) baling d), for various ranges of the variable x. We begin with the result 

'fcaiin g (* 5 d) = £ e-*«t-V(-P, ~xt) |E e~ iM - ( 47 rf)- d/2 - 1 J dt 

= (4^)- d/2 V/ e -^-( d / 2 +f) 7 *(-p,-xt)e- fc2 / 4t ^ (II) 
where use has been made of the Poisson sum formula and the fact that 

POO 

/ e- xt t- p j* (-p, -xt)dt = (12) 
Jo 

The proof of this last statement is relatively straightforward. Making use of the series representation of the function 
7* one obtains 
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e -«*t-P 7 *(_p } - x t)dt = x~ l+p I e"*t- p 7*(-p, — t)di 

) 

^ 00 P 00 4-n—p 



3C 



oo 

-1+p. 



n=0 

OO 



T (-p)l~kJo (n-p)nl dt 



x - 1+p ~ F ( 1 + n _ p ) 
r (-rf ^ (n-p)n! 

(13) 



The last equality in (|^) holds because 



T (-P)h (n-p)n\ (M) 

From Eq. ( pl| ) it is clear that, if a; — > oo and xt — O(l), then the contributions arising from t ~ 1/x will be 
exponentially small in x because of the term going as e~ k / 4t . Therefore, the leading order contributions in the 
regime x — > oo will be generated by the asymptotics of 7*(— p, xt) in the regime xt 3> 1. Making use of the identity 
{§ 

|^ [e*x° 7 * («, a:)] = e x x a - V (« - «, x) (15) 
One can write (1 + xd/dx)I^ caling (x, d) in the following form 

x [(p + 1) 7 * (-p, -xt) + 7*(-p - 1, -xt)] dt (16) 

We now need the asymptotics of the function 7*(— p, —x) for x — ► oo. The leading order behavior of the function in 
this limit follows form the relationship 

7 *(-a,-z)= 1 M(-a,l-a,x) (17) 
T(l - a) 

where M{a, 6, z) is the Kummer function, and the corresponding asymptotic behavior of that function is p6[ 

M(a 1 b,z) = ^-e z z a - b [l + 0{\z\- 1 )], 5iz > 0, \z\ oo (18) 
T(a) 



One then obtains 



7*(-a, -x) ~ — ^— e 1 ! -1 ( 1 + — ) , o<l, x»l (19) 



r(— a) \ x 

where C is a constant. In order to determine C we make use of the identity [Bq 



e 

7*(a,-x) = -x7*(a + l,-x) + — ; — — (110) 

1 (a + 1) 



Making use of this equation and (||) we find that C = a, i.e. 

7*(-o, -x) ~ — !— e 2 ^ 1 f 1 + -) , a < 1, x > 1. (Ill) 
I (—a) V x/ 



\-a) 

Inserting this result into (M), we obtain 



x- 1 (112) 
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where 



l+P 



r(- P ) (4tt) 



/•OO 



+d/2) e - fe V4* & 



(113) 



The above holds when < p < 1. It is easy to show that C p > 0. If the integral in (113) is performed, C p is recast in 
the form 



(i+p)4 1 +"r(i+p + rf/2) 



y 1 

ud+2( P +i) 



T d/2 r(-p) ^ fcd+2( P +i) 1 



In terms of Madelung type constants C(d \a) one can rewrite C p in the form 
or, equivalently in terms of Epstein zcta function Z this constant is 



(i + p)4i+pr(i+p + d/2) 



T d/2 



r(-p) 



(d,-+ P + i). 



(114) 



(115) 



(116) 



For the case d = a = 3 (then p = 1/2) the numerical evaluation gives Ci/ 2 = 10.216. 

When the parameter p is equal to zero, appropriate to the case of short-range interactions, the asymptotic form of 
interest is of the function 



dte~ 



k 

pOO 

(Av)- d/2 V / r^e-^'^dt - l/o 
, Jo 



fc#0 



1 



We now make use of the asymptotic form of the modified Bessel Function 

K v {x) ~ \j^, e ~ X i x > !• 

We immediately find that, for x ^> 1, 



i 



d%/2 



^scaling 

l 



a .(d-3)/4 e -Vi 



(117) 



(118) 



(119) 



Let us now derive the asymptotic form of (1 + xd '/ 'dx)I^ cahne .(x , d) for x -C 1- To that end one needs only to note 
that 



e x l*(a,x) 



r(a + i) 



(120) 



when x — > 0. Then, as a; — ► 0, 

Scaling - Scaling (M) 



r( 



r(i-p)7 



(47Ti)-' i / 2 ^ e " fc2 / 4t -l 



(7f 



(121) 
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The integrands on the right hand side of Eq. (121) are well-defined for p < 1 and 2 < d < 4, both at the lower and 
upper bounds of integration. We will denote the above constant by D p . Then, 



d \ 

1 + 0<p<l, 



x -> 0. 



When p = 0, it is straightforward to show that 

^scaling { x i — ^scaling (0> d) = Do- 
It is easy to show that 



D p = 



(122) 



(123) 



(124) 



Since the Madelung constants C(d\a) are negative for d/2 > a, we obtain that D p < (for p < 1). The numerical 
evaluation for d = a = 3 (then p = 1/2) gives -D1/2 = —0.452 while -Do = —0.226 (which is consistent with 
C(3|l) = -5.029 and C(3|l/2) = -8.914, respectively). 
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FIG. 1. The contour utilized in the contour integration identity (E2) for the summand in the correction to the equation of 
state. 
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FIG. 2. The deformation of the contour in Figure |l]that leads to the new expression (E3) for the summand in the correction 
to the equation of state. 



24 



2 

1 . 75 
1 . 5 
1 .25 
1 

. 75 
. 5 
0.25 




2 4 6 8 10 



FIG. 3. The sum 1/(2: + fc 2 ), as a function of zL 2 , divided by L 2 . The infinite system limit to the sum has been removed. 
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FIG. 4. The distortion of the contour displayed in Figure hi that leads to the expression (F2) for the equation of state 
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